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Effective field theory of the disordered Weyl semimetal 
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In disordered Weyl semimetals, mechanisms of topological origin lead to the protection against 
Anderson localization, and at the same time to different types of transverse electromagnetic response 
- the anomalous Hall, and chiral magnetic effect. We here apply field theory methods to discuss 
the manifestation of these phenomena at length scales which are beyond the scope of diagrammatic 
perturbation theory. Specifically we show how an interplay of symmetry breaking and the chiral 
anomaly leads to a field theory containing two types of topological terms. Generating the uncon¬ 
ventional response coefficients of the system, these terms remain largely unaffected by disorder, i.e. 
information on the chirality of the system remains visible even at large length scales. 
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Weyl semimetals are paradigmatic examples of gap¬ 
less topological condensed matter systems. A Weyl 
semimetal comprises an even number of linearly dis¬ 
persive band touching points embedded in a three- 
dimensional Brillouin zone. The presence of these hot¬ 
spots implies a response to perturbations that is interme¬ 
diate between that of metals and insulators (see Ref. [T] 
for review). This ‘semi-metallicness’ also signifies in the 
physics of the disordered system di]: on the one hand, 
the vanishing of the nodal density of states weakens dis¬ 
order scattering cross sections; on the other hand, suf¬ 
ficiently strong disorder will generate a finite band cen¬ 
ter density of states to eventually overpower the above 
effect. It has been shown |3] that the above competi¬ 
tion manifests itself in the presence of a critical disorder 
strength below/above the system flows towards a clean 
fixed point/a regime of strong impurity scattering. It is 
the purpose of the present paper to derive and discuss 
the effective theory describing the latter phase at length 
scales exceeding the system’s scattering mean free path. 

At large length scales, impurity scattering will render 
the motion of individual excitations diffusive which sug¬ 
gests that the system will end up in the universality class 
of the 3d Anderson metal (i.e. above the phase transition 
point separating a 3d metal from an insulator.) This ex¬ 
pectation is, in fact, a certainty given that a single Weyl 
node may be interpreted as an effective surface theory of 
a bulk 4d topological insulator; finite conduction is pro¬ 
tected by topology. At the same time, topology implies 
a number of differences distinguishing the Weyl system 
from a generic metal: First, an individual Weyl node 
breaks parity symmetry, and it is known (e.g. from the 
example of the d-wave superconducting quasiparticle sys¬ 
tem 0 ) that the breaking of discrete symmetries is gener¬ 
ally remembered, even in the presence of strong disorder. 
Indeed, we will find that the low energy theory of individ¬ 
ual nodes system contains a parity breaking non-abelian 
Chern-Simons (CS) term, which describes the survival of 
the so-called chiral magnetic effect (CME) [7l |8] in the 
disordered environment. Second, it has been shown that 


a system comprising two Weyl nodes separated in mo¬ 
mentum space shows an anomalous Hall effect (AHE) [5] . 
Within the field theoretic framework below, this effect 
will derive from a 3d extension of a two-dimensional topo¬ 
logical 0-term, familiar from the theory of the quantum 
Hall effect. 

Field theory — In the following, we will derive an effec¬ 
tive feld theory describing these structures. Our starting 
point is the bi-nodal Hamiltonian (cf. Fig. 

H = v%(j^ + {v$ + pl)+ V{-x), (1) 

where ^ = k • cr, cr is a vector of Pauli matrices, k the 
vector momentum operator, and v a characteristic ve¬ 
locity. The Pauli matrix cTg acts in a two-component 
space discriminating between two nodes split by a vec¬ 
tor 2b = 26e3 in momentum space and an increment 2/i 
in energy. The model is coupled to disorder by a Gaus¬ 
sian distributed potential P(x) with variance 70. We dis¬ 
criminate between disorder correlated over length scales 
>6“^, which is soft in the sense that the two Weyl nodes 
are not coupled by impurity scattering, and the opposite 
case of short range correlated disorder mixing the nodes. 
Implicit to the model is a high momentum cutoff |k| < A 
limiting the range of linearizability of an underlying lat¬ 
tice model. 

To access the transport properties of the system at 
energies e, we introduce a supersymmetric generating 
functional Z[a\ = J £)(■!/:, ^/>) exp(—iS'[i/j, i/;]), with ac¬ 
tion = —i f d^xif(e + — H)'ip. Here ip = 

is a sixteen-component field of integration 
variables where, n = 1, 2 labels the two nodes, i = 1,2 
the components of a Weyl spinor, s = ± distinguishes 
between advanced and retarded (ar) Green functions, 
i.e. {t3)ss' = s 6 ss' such that integration over ips=± gen¬ 
erates matrix elements of = {e zt i 6 — H)~^, and 
a = 1,2 defines a two-dimensional space (bf-space) of 
complex commuting/Grassmann variables ip^^'^. In the 
commuting sector, ip^ = ip^^, while the commuting vari¬ 
ables ip"^, ip'^ are independent Grassmann variables. From 
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Z^ transport observables may be computed by introduc¬ 
ing suitably defined source fields a, which are presently 
suppressed to keep the notation simple. 

To explore the infiuence of disorder on the system, 
we integrate over V to generate the quartic contribution 
^ f dx The fate of this nonlinearity under changes 

of the cutoff A has been studied 0 [in] by evaluating 
the results in one- and two-loop renormalized perturba¬ 
tion theory in 2 -|- e dimensions El at e = 1. It has 
been found that for bare amplitudes larger than a crit¬ 
ical value 7 * = tt^u^/A the effective disorder strength 
increases under renormalization. To describe the pertur- 
batively inaccessible regime beyond the scattering mean 
free path, A“^ = I ~ we decouple the nonlinear¬ 

ity by a supermatrix field B, and integrate over the then 
quadratic fermions El- As a result, we obtain the effec¬ 
tive action ^[i?] = ^ f d^xstrB^ — strln(G'[i3]), where 

G[B] = (e -I- iSr^ — Hq — B)~^, Hq is the clean Hamil¬ 
tonian and str is the supertrace El- The difference be¬ 
tween the cases of hard and soft disorder, resp., is that in 
the former/latter case the two nodes couple to the same 
{B = H ® P)/independent {B = bdiag(Bi,H 2 )") ma¬ 
trix fields. For definiteness, we first consider the soft 
case, for which the two nodes can be discussed sepa¬ 
rately; the effect of impurity mixing can be described 
by a locking Bi = B 2 at any later stage. Focusing on 
node n = 1 and writing B = Bi for notational simplic¬ 
ity, we proceed by subjecting the action to a mean field 
analysis. A variation of the action yields the equation 
H = ytr G'(x,x; [B]), which is solved [TT] by the diagonal 
ansatz, B = describing the ‘spontaneous symme¬ 

try breaking’ of the infinitesimal (5 a finite impurity self 
energy n = K(e). Specifically, at e = 0 (semimetal) one 
finds K = (2/7r)uA(l — 7 */ 7 ) E]> while far away from 
Weyl node (metal) k = 'jtti/, where ly = e^/27r^u^ is the 
clean density of states. 

Turning to the analysis of fluctuations around the 
mean field, we focus on configurations B = inQ = 
whose action vanishes in the limit of slowly 
varying fluctuation matrices T(x) —>■ const. Here, T G 
G/H takes values in a Goldstone mode coset space where 
G = GL(2|2) is the group of 4 x 4 invertible super¬ 
matrices and H = GL(2|2) the subgroup of matrices 
k commutative with the mean field, [fc,r 3 ] = 0. Since 
str(Q^) = 0, the action we now need to consider reads 
5 'o[Q] = stTln{e — + inQ) . For later reference, we note 

that the action has two continuous symmetries: a global 
symmetry Q —>■ TqQT^^, where To S G is constant, and 
a local gauge symmetry T Tk, k = k{x.) G H. 
Effective field theory — Our goal is to expand the ac¬ 
tion in gradients T~^diT = Ai. One might be tempted 
to start this program with a similarity transformation, 
'S'o[Q] = strln(T“^(e—u|f-|-i«:Q)T) = str ln(e— 
iKTa) = S'[A], which may then be expanded in powers of 
the ‘non-ablian gauge field’ A. However, due to the noto- 
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Figure 1. (Color online) Schematic of two Weyl nodes split 
in energy and momentum by 2y, and 2b, resp. 

rious anomaly, this operation is invalid, the action needs 
to be regularized first. Following a strategy previously 
applied to the 2 d d-wave superconductor, we regularize 
by subtraction of a term S'[(5] = S'o[(5] — where Sri 

differs from Sq by a replacement k ^ r] \i 0 and setting 
e = 0. In the limit 77 -i- 0 the Q-dependence of Sri[Q] 
drops out so that = 0. On the other hand, for 

large momenta v\k\ ^ q and fixed rj, the two action con¬ 
tributions cancel against each other, i.e. 5'[Q] is UV reg¬ 
ularized (in a gauge-invariant way.) The similarity trans¬ 
formation may now safely be applied to both S'o,?;, so that 
we obtain an effective action S'[Q] = (5'o[A] — 
where the subscript ‘reg’ means that only UV finite con¬ 
tributions to the subsequent expansion in A are to be 
kept. Notice that in the language of the A-fields, the 
local gauge invariance is no longer manifest. Rather, 
Ai -G k~^{A -f kdik~^)k transforms as a non abelian 
gauge field, and invariance of the action becomes a non¬ 
trivial consistency condition. 

In the expansion of the action, we keep terms of order 
two (G(9A, A^)) and three (G(A^,9A)) derivatives. To 
second order we obtain the result 

*5d[A]=-^^ I d^XStff[Ar,Tsf), 

<5'top[A] = J d^x str{rsdiAj), (2) 

where the longitudinal and Hall conductivity of node 1, 
and al,y are determined by the microscopic model 
parameters as discussed below. We note that the ac¬ 
tion ([^ affords the manifestly gauge invariant refor¬ 
mulation S'd[Q] = — ^ f d^xstr(dQ^) and S'top[Q] = 
J dx str{QdiQdjQ). The first of these contribu¬ 
tions has been constructed in Ref. El on phenomeno¬ 
logical grounds within a non-regularized framework. (In 
our current approach regularization plays a vital role; the 
‘cross term’ AT 3 AT 3 of the commutator [A, t^]'^ is weighed 
with an UV divergent ‘fermion bubble’. Only after reg¬ 
ularization it combines with the finite coefficient of A^ 
to a gauge invariant commutator.) Before discussing the 
physics of these expressions, we complete the derivation 
of the action and consider terms of cubic order in A. 
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The terms of O(A^) are the ‘triangle graphs’ pervasive 
in the theory of (2 + 1) or (3 + 0) dimensional relativistic 
gange theories. On general grounds [14] we expect the 
appearance of a Chern-Simons action at this order. A 
straightforward if lengthy calculation indeed yields the 
result 

5cs[^] = 56s[^] + *56i[A], (3) 

jf-ijk r 

^cs[^] = ^ E d?xsiY{A,P^d,AuP^), 

jfijk r 

S=± 

where P^ is a projector on advanced/retarded indices. 
Apart from the presence of these projector matrices, this 
has the characteristic structure of a non-abelian CS term. 
(However, ^ does not define a ‘real’ CS action, inasmuch 
as A does not describe a field-gauge coupling, but repre¬ 
sents the nonlinear u-model field itself; the situation is 
conceptually similar to that considered in Refs. (ISIIISl.) 

The CS action does not afford a representation in 
terms of Q-fields, which reflects the lack of com¬ 
plete gauge invariance of this action piece M- it is 
straightforward to verify that under a gauge transfor¬ 
mation by fc = bdiag(A:_|_, A:_S H, and for a ficti¬ 
tious system without boundaries, the CS action trans¬ 
forms as Scs[A] + •S'toplfc], where Stop[k] = 

^ s f d^xstr(k~^dks)^^. The integral yields a 

quantized value, viz. 247r^ x Ug, where rig is the is wind¬ 
ing number of a configuration kf G SU(2) in three di¬ 
mensional space. (The non-compact bb sector of the su¬ 
permatrices k is topologically empty.) For ‘large’ gauge 
transformations with non-vanishing winding numbers, 
the CS action changes by a factor f7r(n+ -|- n_). The 
origin of this phenomenon was explained in the clas¬ 
sic reference El. where it was shown that these factors 
get canceled by a gauge anomaly of the regulator action 
Sfj ^ Sjj + i7r(n+ -I- n_), where the appearance of the 
extra terms is caused by zero-crossings of the regulariz¬ 
ing Dirac operator under a large gauge transformation. 
Overall the action is gauge invariant. 

Discussion — The action S'^-I-S'top of the diffusion modes 
A appeared for the first time in connection with the 
multi-layer quantum Hall effect m, a system conceptu¬ 
ally similar to the present one if the 3-direction of Weyl 
node splitting is interpreted as the stacking direction of 
a layered system of 2d class A topological insulator lay¬ 
ers |5]. The action controls the fluctuations of diffu¬ 
sion modes in terms of the dimensionless coupling con¬ 
stant g^x = 'Xxx^~^ Within the framework of our gra¬ 
dient expansion we find Uxx = (e^ + 3/c^)/67rKn, which 
simplifies to = K/27rn at the Weyl node and asymp¬ 
totes to the Drude conductivity = n ^/37 at higher 
energies e k. At the nodes, e = 0, and at bare length 
scales A ~ characteristic for the ballistic/diffusive 


crossover, the conductance gxx takes values of 0 ( 1 ), close 
to but larger m than the critical value g* marking the 
3d Anderson transition. 

The bare coefficient a^y multiplying the action S'top 
is the contribution of node 1 to the Hall conductiv¬ 
ity of the system at crossover length scales to the dif¬ 
fusive regime. As in the derivation of the field the¬ 
ory of the quantum Hall effect m we obtain axy as 
a thermodynamic coefficient axy = V~^dBN, where 
V = LzL^ is the volume of the sample, and the deriva¬ 
tive probes the sensitivity of the number of states N = 
— fieImtr(C?+(e)) below zero energy to the pres¬ 

ence of an external magnetic field in 3-direction. Within 
the framework of the gradient expansion, G~*' is the re¬ 
tarded Weyl node Green function coupled to disorder 
in the mean field/self consistent Born approximation. 
Following the strategy of Ref. [3], we evaluate this ex¬ 
pression by interpreting the quantized values of the 3- 
momentum, m„ = v(b + k^^n), ks = 2’nnL~^ as masses 
entering effective 2d layered Dirac Hamiltonians govern¬ 
ing the planes perpendicular to the 3-direction. Within 
this interpretation, axy = {‘2-'x/Lz)'^^axy,n is obtained 
by summing over the contribution of the layers, where 
'Xxy,n = Cn!2 t: is given by the disorder averaged Chern 
number, C„ = (l/27r) X]o'=± arctan [(to„ -|- ae)/^] -I- 1/2, 
of the nth layer. Since only changes of at zero cross¬ 
ings of the effective masses can be unambiguously de¬ 
termined from the linearized theory, we fix the absolute 
value of the sum by the condition axy = 0 at & = 0, 
which follows from matching to the band structure of 
the bi-nodal Weyl system. In the limit L^b ^ 1, this 
gives axy = &/27r for both nodes, irrespective of the en¬ 
ergy e or the disorder strength k. Both for soft and hard 
disorder, the two contributions to the Hall conductivity 
add, and we obtain axy = b/n, a result known as the 
anomalous Hall effect (AHE). 

The CS contribution to the action accounts for the 
thermodynamic response of the system to imbalances 
(b, g) between the nodes. That this term probes equi¬ 
librium properties follows from the presence of the pro¬ 
jector matrices P® which prevent coupling between the 
retarded and advanced sector of the theory (a necessary 
ingredient to any type of dynamic response.) A finite 
equilibrium response is obtained if we couple the sys¬ 
tem to an external field a = {o^}, where ai = ^e^ijX^, 
i = 1,2 represents an external magnetic field Be^, and 
03 = a(x)a^^ (g) T 3 is a source field. The latter is de¬ 
fined in such a way that differentiation of the parti¬ 
tion function, -^5a(^)Z[a] = -^Im (tr[G+(x,x)CT3]) = 
{j 3 {x.)S{€ — H)) = yields the contribution of states at 
energy e to the equilibrium value of the 3-current density 
of node 1. We compute this expression by adding the 
external field to the internal one, A ^ A + a, and substi¬ 
tuting this configuration into the CS action. In the sim¬ 
plest approximation A = 0 (for the above ‘equilibrium’ 
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choice of source terms fluctuation corrections around the 
A = 0 vanish by supersymmetry anyway), we then ob¬ 
tain 5'cs[a] = —^ / dxa{x), and hence To 

obtain the full response of the system, we need to add 
the (opposite) contribution of the second node and inte¬ 
grate over filled energy states up to some Fermi energy 
ep. Taking into account that the existence of a bare lin¬ 
earization cutoff A implies a cutoff |e| < |eo ±/i|, eo = vA 
for the accessible energy states (cf. Fig.[^, this leads to 
J 3 = dej 3 ,e - dej 3 ,e = i.e. an equilib¬ 

rium current proportional to an external magnetic held, 
the so-called chiral magnetic effect (CME) [3 [8] (for a 
discussion how the nonvanishing of this expression may 
be understood from the perspective of the Fermi-liquid 
theory, see Ref. [10]). 

Renormalization - What happens if short distance fluc¬ 
tuations in the held theory are integrated out to probe 
the physics at length scales beyond the ballistic/diffusive 
crossover regime? An answer to this question has been 
formulated in Ref. HU within the framework of two loop 
renormalized perturbation theory for the dimensionless 
coupling constants = a^^A~^ of the model. The re¬ 
sult 


dgxx 1 

d In I ^'x'^gxx ’ 


dgxy _ 

di^i = 


(4) 


states that the longitudinal conductance scales accord¬ 
ing to the predictions of one-parameter scaling theory 
(unaffected by the Hall conductance) towards Ohmic be¬ 
havior gxx 1. The Hall conductance shows linear 

scaling, gxy oc I, which means that the AHE remains 
unrenormalized by disorder cr^y = const, even at large 
length scales (in contrast to the Hall conductivity of a 
two-dimensional system which is renormalized by instan- 
ton fluctuations m)- Finally, the coupling constant of 
the CS action is fixed by gauge invariance, and fluctua¬ 
tion corrections to (ja) b vanish by supersymmetry. This 
means that within the framework of our theory the CME, 
too, fully protected against renormalization by disorder. 

The disorder insensitivity of the topological response 
coeffecients holds regardless of whether we are probing 
the semimetallic Weyl nodes e, /i ^ k, or the metallic 
physics far away from them |e ± /r| ^ k. The essen¬ 
tial difference between the two situations lies in the bare 
and renormalized values of the longitudinal conductance 
gx:x- in the former/latter case, gxx is initially small/large 
to begin with. However, in either case, gxx increases, 
and asymptotes to Ohmic behavior at large length scales. 
While the system then behaves similar to a three di¬ 
mensional metal, the preserved non-vanishing of its two 
transverse transport coefficients betrays the underlying 
presence of two Dirac nodes. 

Summarizing, we have microscopically derived a su¬ 


persymmetric held theory description of disordered Weyl 
semimetals and metals at length scales exceeding the 
mean free path. The profile of the theory is essentially de¬ 
termined by an interplay of symmetry conditions and the 
chiral anomaly. This mechanism stabilizes metallic be¬ 
havior at large length scales, along with various disorder- 
insensitive response coefficients of topological origin. 
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